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It is known that in many cases distributions of exponential integrals of Levy processes are in- 
finitely divisible and in some cases they are also selfdecomposable. In this paper, we give some 
sufficient conditions under which distributions of exponential integrals are not only selfdecom- 
posable but furthermore are generalized gamma convolution. We also study exponential integrals 
of more general independent increment processes. Several examples are given for illustration. 
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1. Introduction 

Let (£,Tj) = {(^t,Vt)it > 0} be a bivariate cadlag independent increment process. In most 
cases, (£, n) is assumed as a bivariate Levy process, but we also treat more general cases 
where £ or r\ is a compound sum process, which is not necessarily a Levy process but is 
another typical independent increment process. Our concern in this paper is to examine 
distributional properties of the exponential integral 



provided that this integral converges almost surely. More precisely, we are interested in 
when C(V), the law of V, is selfdecomposable and moreover is a generalized gamma 
convolution. 
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We say that a probability distribution p, on K (resp. an R-valued random variable X) 
is sclfdecomposablc, if for any b > 1 there exists a probability distribution /it, (resp. a 
random variable Y& independent of X) such that 

fi = D b -i (fi) * fx b (resp. X = b~ l X + Y b ), 

where D a (p) means the distribution induced by D a (p)(aB) := p(B) for B e Q3(R), * is 

the convolution operator and = denotes equality in distribution. Every sclfdccomposable 
distribution is infinitely divisible. Some well-known distributional properties of nontrivial 
selfdecomposable distributions are absolute continuity and unimodality (see Sato [15], 
pages 181 and 404). 

First, we review existing results on C(V). Bertoin et al. [3] (in the case when r\ = {774} 
is a one-dimensional Levy process) and Kondo et al. [9] (in the case when r\ is a multi- 
dimensional Levy process) showed that if £ = {£ t } is a spectrally negative Levy process 
satisfying limt->.oo£t = +00 a.s. and if the integral (1.1) converges a.s., or equivalently, if 
/ R log + |y|^(dy) < 00 for the Levy measure v n of rji, then C(V) is selfdecomposable. 

On the other hand, there is an example of noninfmitely divisible C(V), which is due 
to Samorodnitsky (see Kliippelberg et al. [8]). In fact if (£t,Vt) = (St +at,t), where {S t } 
is a subordinator and a > some constant, then the support of C(V) is bounded so that 
C(V) is not infinitely divisible. 

Recently, Lindner and Sato [11] considered the exponential integral 

[ exp(-(logc)iV t _)dy t = f C Nt -dY t , c>0, 

where {(N t ,Y t )} is a bivariate compound Poisson process whose Levy measure is concen- 
trated on (1,0), (0, 1) and (1, 1), and showed a necessary and sufficient condition for the 
infinite divisibility of C(V). They also pointed out that C(V) is always c _1 -decomposable, 
namely there exists a probability distribution p such that [i = D c -i(p.) * p. Note that a 
c _1 -decomposable distribution is not necessarily infinitely divisible, unless p is infinitely 
divisible. In their second paper (Lindner and Sato [12]), they also gave a condition under 
which C(V), generated by a bivariate compound Poisson process {(7V t ,y t )} whose Levy 
measure is concentrated on (1,0), (0,1) and (l,c _1 ), is infinitely divisible. 

For other distributional properties of exponential integrals, like the tail behavior, see, 
e.g., Maulik and Zwart [13], Rivero [14] and Behme [2]. 

In this paper, we focus on "Generalized Gamma Convolutions" (GGCs, for short) to 
get more explicit distributional informations of V than selfdecomposability. 

Throughout this paper, we say that for r > and A > a random variable "f r .\ has a 
gamma(r, A) distribution if its probability density function / on (0,oo) is 

f( X ) = -?- X r-l e -^. 

r(r) 

A gamma(l, A) distribution is an exponential distribution with parameter A > 0. When 
we do not have to emphasize the parameters (r, A) , we just write 7 for a gamma random 
variable. 
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The class of GGCs is defined to be the smallest class of distributions on the positive 
half line that contains all gamma distributions and is closed under convolution and weak 
convergence. By including gamma distributions on the negative real axis, we obtain the 
class of distributions on M which will be called "Extended Generalized Gamma Convolu- 
tions" (EGGCs, for short). We refer to Bondesson [4] and Steutel and van Harn [16] for 
many properties of GGCs and EGGCs with relations among other subclasses of infinitely 
divisible distributions. 

One well-known concrete example of exponential integrals is the following. When 
= (Bt + at,t) with a standard Brownian motion {B t } and a drift a > 0, the law 
of (1.1) equals £(1/(27)) which is GGC (and thus is selfdecomposablc). 

When choosing £ to be deterministic, that is, (£,t,Vt) = (tiVt), the exponential integral 
(1.1) is defined and is an EGGC if and only if 77 admits a finite log-moment (needed 
for the convergence) and C(rji) is included in the Goldie-Steutel-Bondesson Class, a 
superclass of EGGC as defined for example, in Barndorff-Nielsen et al. [1]. This fact 
follows directly from [1], equation (2.28). In the same paper, the authors characterized the 
class of GGCs by using stochastic integrals with respect to Levy processes as follows. Let 
e{x) = J x u~ 1 e~ u du, x > 0, and let e*(-) be its inverse function. Then £(/( ^ e*(s) dr] s ) 
is GGC if 77 is a Levy process and £(771) has a finite log-moment. 

In this paper, via concrete examples, we investigate distributional properties of expo- 
nential integrals connected with GGCs. 

The paper is organized as follows. In Section 2, we give some preliminaries. In Sec- 
tion 3, we consider exponential integrals for two independent Levy processes £ and 77 such 
that £ or 77 is a compound Poisson process, and construct concrete examples related to 
our question. In the special case that both £ and 77 are compound Poisson processes, we 
also allow dependence between the two components of (£,77). In Section 4, we consider 
exponential integrals for independent increment processes such that £ and 77 are inde- 
pendent and one is a compound sum process (which is not necessarily a Levy process) 
while the other is a Levy process. 

2. Preliminaries 

The class of all infinitely divisible distributions on K (resp. R+) is denoted by (resp. 
7(K+)). We denote the class of selfdecomposable distributions on K (resp. K+) by L(K) 
(resp. L(R+)). The class of EGGCs on M (resp. GGCs on K+) is denoted by T(R) (resp. 
T(R + )). The moment generating function of a random variable X and of a distribution 
H are written as Lx and L M , respectively. If X is positive and /i has support in R + , Lx 
and L M coincide with the Laplace transforms. 

We are especially interested in distributions on R + . The class T(R+) is characterized 
by the Laplace transform as follows: A probability distribution /j, is GGC if and only if 
there exist a > and a measure U satisfying 




and 
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such that the Laplace transform L M (z) can be uniquely represented as 

hJu)= [ e- ux f i(dx)=expl-au+ [ logf ) U(dx)\. (2.1) 

-'[Coo) I J (0,oo) \X + Uj J 

Another class of distributions which we arc interested in is the class of distributions 
on R + whose densities are hyperbolically completely monotone (HCM, for short). Here 
we say that a function f(x) on (0,oo) with values in R + is HCM if for every u > 0, the 
mapping f{uv) ■ f(u/v) is completely monotone with respect to the variable w = v + v^ 1 , 
v>0. Examples of HCM functions are x f) {f3 G R), c~ cx (c > 0) and (l + cx)~~ a (c>0,a> 
0). The class of all distributions on R + whose probability densities are HCM is denoted 
by iJ(R+). Note that H(R + ) c T(R+) c L(R + ) C /(M+). 

For illustration, we give some examples. Log-normal distributions are in H(R + ) [4], 
Example 5.2.1. So these are also GGCs. Positive strictly stable distributions with Laplace 
transform L(w) = cxp(— u a ) for a G {1/2, 1/3, . . .} arc in H (R+) [4], Example 5.6.2, while 
they are GGCs for all a G (0, 1] [16], Proposition 5.7. Let Y = exp(7 r .A) — 1. If r > 1, then 
£{Y) is in F(R+), but if r < 1, C(Y) is not in H(R + ) [4], page 88. But C(Y) or equiva- 
lently £(exp(7 rj A)) is always in T(R + ), independent of the value of r [4], Theorem 6.2.3. 
Remark that by treating H(R + ) we cannot replace exp(-f Tt \) — 1 by exp(j rt \). Namely, set 
r = 1 and observe that the probability density function X(x + 1)^ x ^ 1 1[q :OQ )(x) is HCM, 
but the probability density function Ax~ A ~ 1 l[ 1 . oc )(x) is not HCM. It follows from this 
that £(cxp(7i jA ) - 1) is in H(R + ) but £(exp(7i jA )) is not in H(R+). 

In addition, we also investigate the modified HCM class denoted by H(M), which 
gives some interesting examples of C(V) on R. The class H(M) is characterized to be 
the class of distributions of random variables \J~XZ, where C(X) £ H(M.+ ) and Z is a 
standard normal random variable independent of X (see Bondcsson [4], page 115). By 
the definition, any distribution in H(M) is a type G distribution, which is the distribution 
of the variance mixture of a standard normal random variable. Note that H(M) C T(R). 
As will be seen in Proposition 2.1, there are nice relations between H(M) and T(R) in 
common with those of H (R+) and T(R+). 

Here we state some known facts that we will use later. 

Proposition 2.1 (Bondesson [4] and Steutel and van Harn [16]). 

(1) A continuous junction L(w), u > 0, with L(0+) = 1 is HCM if and only if it is the 
Laplace transform of a GGC. 

(2) If C(X) 6 H(R+), C(Y) G T(R+) and X and Y are independent, then C(XY) G 
T(R+). 

(3) Suppose that C{X) G H(M.+), C{Y) G T(R) and that X and Y are independent. If 
C(Y) is symmetric, then C(-/XY) G T(R). 

(4) Suppose that C(X) G H (R) and C(Y) G T(R) and that X and Y are independent. 
If C{Y) is symmetric, then C{XY) G T(R). 

(5) If C(X) G H(R), then C{\X\*) G H(R+) for all \q\ > 2,q G R. Furthermore, 
£(|X|9 sign(X)) G i?(R) /or all q G N, g ^ 2, but not always for q = 2. 
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Remark 2.2. Notice that the distribution of a sum of independent random variables 
with distributions in H(M.+) does not necessarily belong to if(R+). See Bondesson [4], 
page 101. 

Some distributional properties of GGCs are stated in the following proposition [4], 
Theorems 4.1.1. and 4.1.3. 

Proposition 2.3. 

(1) The probability density function of a GGC without Gaussian part satisfying < 
/(o oo) U(du) = f3 < oo with the measure U as in (2.1) admits the representation x@~ 1 h(x), 
where h{x) is some completely monotone function. 

(2) Let f be the probability density of a GGC distribution without Gaussian part sat- 
isfying 1 < J/ Q v U(du) — P < oo. Let k be a nonnegative integer such that k < (3 — 1. 
Then the density f is continuously differentiable any times on (0,oo) 7 and at at least 
k times differentiable with /^(0) = for j < k. 

Examples of GGCs and the explicit calculation of their Levy measure are found in 
Bondesson [4] and James et al. [7]. 

Necessary and sufficient conditions for the convergence of (1.1) for bivariate Levy 
processes were given by Erickson and Mailer [5]. More precisely, in their Theorem 2.1, 
they showed that V converges a.s. if and only if for some e > such that A^(x) > for 
all x > e it holds 



Here A^(x) = <Zf + z^((l, oo)) + , ^((y, oo)) dy while (Y<x,vx,ax) denotes the Lcvy - 
Khintchinc triplet of a Levy process X . 

3. Exponential integrals for compound Poisson 
processes 

In this section, we study exponential integrals of the form (1.1), where either £ or i] is a 
compound Poisson process and the other is an arbitrary Levy process. First, we assume 
the two processes to be independent, later we also investigate the case that (£,77) is a 
bivariate compound Poisson process. 

3.1. Independent component case 

We start with a general lemma which gives a sufficient condition for distributions of 
perpetuities to be GGCs. 



lim £ t = 00 



a.s. and 




(2.2) 
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Lemma 3.1. Suppose A and B are two independent random variables such that C(A) £ 
H(R+) and C(B) £ T(R + ). Let (Aj, Bj),j = 0, 1, 2, . . . , be i.i.d. copies of (A,B). Then, 
given its a.s. convergence, the distribution of the perpetuity Z := X^odliLo 1 Ai)Bk be- 
longs to T(M_|_). Furthermore, if C{A) £ H(R), C(B) £ T(R) and C{B) is symmetric, 
then C{Z) £ T(R). 



Proof. If we put 

n /fc-1 \ 

then we can rewrite 



Z n = B + A {B 1 +A 1 (B 2 + --- + A„_ 2 (5 n _! + VA) • ■ ■))• 

Since A n ~i and B n are independent, e H(R + ) and C(B n ) £ T(R+), we get 

£(A n _i£? n ) S T(R+) by Proposition 2.1. Further and A n ^\B n are independent 

and both distributions belong to T(R+) and hence L(B n -\ + A n _\B n ) £ T(R+). By 
induction, we can conclude that C(Z n ) £ T(R + ). Since the class T(R+) is closed under 
weak convergence, the first part follows immediately. A similar argument with (4) in 
Proposition 2.1 gives the second part. □ 



Case 1: The process £ is a compound Poisson process 

Proposition 3.2. Suppose that the processes £ and r\ are independent Levy processes 
and further that £ t = * s a compound Poisson process with i.i.d. jump heights 

Xi,i = 1,2, ... , such that < -E[A"i] < oo, C(e~ Xl ) £ i/(R+) and r\ has finite log-moment 
£'log + 1 771 1 and it holds C(r) T ) £ T(R+) for an exponential random variable r independent 
of rj having the same distribution as the waiting times of N . Then the integral (1.1) 
converges a.s. and 

C[ [ e-^-drjt) GT(R+). 

\J(P,oo) J 

Furthermore, if £(c~ Xl ) £ H(M), C(r/ T ) £ T(R) and C{r) T ) is symmetric, then C(V) £ 
T(R). 

Proof. Convergence of the integral follows from (2.2). 

Set T = and let Tj,j = 1, 2, . . . , be the time of the jth jump of {N t ,t > 0}. Then we 
can write 

p OO p OO p 

/ e -« f -d77 t =V/ e-S=i x 'dr, t = yc-S=i x ' / dn t 

j=0 \i=l / 
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where £° =1 = 0, JlLi = !> ^ = e x * and = J (T . iT . +l] d?y t = ^t j+1 -t- n ow 
Lemma 3.1 yields the conclusion. □ 

In the following, we first give some examples for possible choices of £ fulfilling the 
assumptions of Proposition 3.2 and then continue with examples for r\. Hence, any com- 
bination of them yields an exponential integral which is a GGC. 

Example 3.3 (The case when X\ is a normal random variable with positive 
mean). We see that £(c~ Xl ) is log-normal and hence is in H(R + ). 

Example 3-4 (The case when X\ is the logarithm of the power of a gamma 
random variable) . Let Y = jr.x and X\ = clogY" for cgl. Recall that £(logF) G 
T(R) and so C(cX t ) = C(logY c ) G T(R) for cel. Note that E[X{\ = ap(X), where 
ip{x) the derivative of logL(x). If we take c G R such that cifj(\) > 0, we conclude that 

Example 3.5 (The case when X\ is logarithm of some positive strictly stable 
random variable) . Let X\ = logF be a random variable, where Y is a positive stable 
random variable with parameter < a < 1. Then X\ is in the class of EGGCs when 
a = 1/n, n = 2, 3, . . . (see Bondesson [4], Example 7.2.5) and 

a uXn _ z^u! _ r(l-u/a) 



E[e uXl ] = E[Y u ] = 



T(l-u) 



It follows that E[X X ] = -^V(l) + ^(1) = (1 - l/a)^(l) > and £( C - Xl ) = Z^F" 1 ) G 
H{R + ) by [4], Example 5.6.2. 

Example 3.6 (The case when X\ is the logarithm of the ratio of two expo- 
nential random variables) . Let X\ = log(Yi/Y2), where Yj,j = 1,2, are independent 
exponential random variables with parameters Xj > 0, j = 1,2. The density function of 
Xi is given by [4], Example 7.2.4, 

/(x)= S(a!,A 2 )(1 + c-)^' x ^ 

where B(-, ■) denotes the Beta-function. Now if 2£[Xi] > we can set X\ to be a jump 
distribution of £. It is easy to see that £(e~ Xl ) = C(Y 2 /Y 1 ) G H(R + ) since £(Yj) G 
ff(M+). 

Example 3.7 (The case when r\ is nonrandom). When rjt = t, it holds that C(rj T ) = 
£(r)GT(R + ). 

Example 3.8 (The case when n is a stable subordinator) . The Laplace transform 
of B := r\ T with r = 71 a is given by (see, e.g., Steutel and van Harn [16], page 10) 

Lfl(u) 



A-logL^(u) 
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Now consider rj to be a stable subordinator without drift. Then the Laplace transform 
of rji is given by L 7)1 (u) = cxp{— Therefore, we have 

L B («) = 



A + u Q 

This function is HCM, since i s HCM by the definition and due to the fact that the 
composition of an HCM function and x a , \a\ < 1, is also HCM (see [4], page 68). Thus, 
the Laplace transform of B is HCM by Proposition 2.1 and we conclude that C(rj T ) is 
GGC. 

Remark that if rj admits an additional drift term, the distribution C(B) is not GGC. 
This result was pointed out by Kozubowski [10]. 

Example 3.9 (The case when rj is an inverse Gaussian Levy process). We 

suppose rj to be an inverse Gaussian subordinator with parameters /3 > and 8 > 0. The 
Laplace transform of rjt is 



L„ t («) = atp(-6t(y/P 2 + 2u- 0)). 
Now by choosing the parameters satisfying A > 8/3, we have, for B = rj T , 

h B (u) : 



A - 8/3 + S^/3 2 + 2u 

This function is HCM by argumentation as in Example 3.8 with a = 1/2 and using 
Property (xi) in [4], page 68. 

Remark 3.10. Although the Levy measure of rj T is known explicitly, it is an open 
question whether the parameter of the exponentially distributed random variable r has an 
influence on the GGC-property of rj T , or not. Examples 3.7 and 3.8 lead to the conjecture 
that there is no influence. So far, no counterexamples to this conjecture are known to 
the authors. 



Case 2: The process rj is a compound Poisson process 

In the following, we assume the integrator 77 to be a compound Poisson process, while £ 
is an arbitrary Levy process, independent of n. We can argue similarly as above to obtain 
the following result. 

Proposition 3.11. Let £ and rj be independent and assume rjt — ^° ^ e a com- 

pound Poisson process with i.i.d. jump heights Yi,i = 1,2, ... . Suppose that E[£i] > 0, 
E\og + 1 771 1 < oo, C{Y\) <= T(R_|_) and £(c~^ T ) e H(R + ) for an exponentially distributed 
random variable r independent of £ having the same distribution as the waiting times 
of N . Then the integral (1.1) converges a.s. and it holds that 

c( f e~t*~ di]t) eT(R + ). 
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Furthermore, if £(e _f -) € H(R), £(Yi) € T(R) and £(Yi) is symmetric, then C(V) € 
T(R). 

Proof. Convergence of the integral is guaranteed by (2.2). Now set To = and let 7}, j = 
1,2,..., be the jump times of {N t ,t > 0}. Then we have 

( . oo oo 



oo / J 



3=1 \i=l / j'=l \i=l 



where Aj = e - ^ - ^-!-' = e~^ T »~ T i-i and = Yj. Remark that the proof of Lemma 3.1 
remains valid even if the summation starts from j = 1. Hence, the assertion follows from 
Lemma 3.1. □ 



3.2. Dependent component case 

In this subsection, we generalize a model of Lindner and Sato [11] and study to which 
class C(V) belongs. 

Let < p < 1. Suppose that (£,77) is a bivariatc compound Poisson process with pa- 
rameter A > and normalized Levy measure 

v{dx,dy) = p6 (dx)p (dy) + (1 - p)S 1 (dx)p 1 (dy), 

where po and p\ are probability measures on (0, 00) and [0, 00), respectively, such that 

/ logydpo(y) < 00 and / logy dpi (y) < 00. 

J (l,oo) J(X, oo) 

For the bivariate compound Poisson process (£,f?), we have the following representation 
(see Sato [15], page 18): 

N t f N t N t \ 

k=0 \fe=0 k=0 ) 

where S^ 1 = = and {Sk = (S^ , S^)}'^ =1 is a sequence of two-dimensional i.i.d. 
random variables. It implies that the projections of the compound Poisson process on R 2 
are also compound Poisson processes. Precisely in the given model, since P(5p' 1 = 0) = p 
and P(Si = 1) = 1 — p, the marginal process £ is a Poisson process with parameter 
(1 - p)X > 0. Note that and S<f ] may be dependent for any k GN. In this case, 
Pi(B) is equal to P(s£ 2 } € B\S[ 1} = i) for i = 0, 1 and B G B(R). 
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Example 3.12. In Lindner and Sato [11], the authors considered the bivariate com- 
pound Poisson process with parameter u + v + w, u, v, w > and normalized Levy mea- 
sure 

v(dx,dy) = — — So(dx)S 1 (dy) + U + W S 1 (dx)( —^—S (dy) + — ^— <5i(dy) 

U + V + W U + V +w \u + w U + W 

So in their setting p = p a = Si and Pi = ^S + ^8 1 . 

In the following theorem, we give a sufficient condition for C(V), given by (1.1) with 
(£,??) as described above, to be GGC. 

Theorem 3.13. If the function ^I^ff is HCM, then C(V) is GGC. 

Proof. Define and M to be the first jump time of the Poisson process £ and the num- 
ber of the jumps of the bivariate compound Poisson process (£, rj) before Tj, respectively. 
Due to the strong Markov property of the Levy process (£,??), we have 



/ exp(-£ 5 _) d?7 s = / exp(-£ s -)dr) s + exp(-£ s _) dr) s 

J(0,oo) J(0,T(] J(T t ,oo] 

Vt ( + cxp(-^ +s -)dr/ s+ T £ 

J(0,oo) 

T) T + exp(-£ T ) / exp(-(£ T +s _ - £r s )) d((r/ (s+T ) - rj ) + 77 ) 

J(0,oo) 5 5 

Vt+^I exp(-f s _)dry s , 

ifO.oo) 



where the process (£, rj) is independent of {(£t, f]t), t < T^} and has the same law as (£, rf). 
Therefore, we have 

L M (u)=L A( (e- 1 ix)L p (u), (3.1) 

with p = C{V) and p denoting the distribution of rj T ^ . Thus, p is e _1 -decomposable and 
it follows that 

oo 

L M («)= l[L p (e- n u). 

n=0 

In the given setting, we have 



' M \ 
<k=0 ) 
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where M is geometrically distributed with parameter p, namely. 

P(M = k) = (1 - p)p k for any k E N Q . 

Hence, we obtain 

L p (w) = E[cxp(-U77 T? )] 

r r / / m 

c(2) 
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E 



E 



exp 



Vfc=0 



7V 



= E[(L P0 ( U )) M L P1 ( U ) 



(1 -p)L pl (u)^(pL po (ii)) 



fc=0 



fe = (1 -p)E pl (u) 
1 - pL Po (u) 



The class of HCM functions is closed under scale transformation, multiplication and 
limit. Therefore, L M (u) is HCM if h p {u) is HCM and hence p, is GGC if p is GGC by 

Proposition 2.1(6). As a result, p is GGC if ^I^ff is HCM. □ 

A distribution with the Laplace transform iSp^ P \ u ) is called a compound geometric 
distribution. It is compound Poisson, because every geometric distribution is compound 
Poisson with Levy measure given by 



u p ({k}) = - 



1 



log(l-p) k + V 



fc = l,2,. 



(see page 147 in Steutel and van Harn [16]). Since HCM functions are closed under 
multiplication, we have the following. 



Corollary 3.14. If pi and the compound geometric distribution of po are GGCs, then 
so is C(V). 



In addition, we observe the following. 



Corollary 3.15. 

(1) For any c > 1, the distribution p c = £(J/ ^ c~^ s ~ d?7 s ) is c _1 -selfdecomposable. 
Thus, in the nondegenerate case it is absolutely continuous or continuous singular 
(Wolfe [17]) and Theorem 3.13 holds true also for p c instead of p. 

(2) If p\ is infinitely divisible, then p c is also infinitely divisible. 

(3) Let _B(R_|_) be the Goldie-Steutel-Bondesson class, which is the smallest class that 
contains all mixtures of exponential distributions and is closed under convolution and 
weak convergence. 

If j~g7^^ is the Laplace transform of a distribution in B(M.+ ), then p c is in B(M + ). 
Moreover, p c will be a c _1 -semi-selfdecomposable distribution. 
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About the definition and basic properties of semi-selfdecomposable distributions, 
see [15]. The proof of (1) is obvious. For (2), remark that a distribution with Laplace 
transform rii! \ as compound Poisson distribution is always infinitely divisible. 
Hence, only p\ has influence on that property. The proof of (3) follows from the charac- 
terization of the class £?(R+) in Chapter 9 of [4] and our proof of Theorem 3.13. 

Example 3.16. Let p x be a GGC, that is, L pi (u) is HCM. Then if jz^Tu) is HCM ' 
p is found to be GGC. For example, if po is an exponential random variable with density 
fix) = bc~ bx ,b > 0, then — / , ■. is HCM. To see this, for u > 0,v > 0, write 

■> \ I ' ' 1— pL P0 (u) ' ' 

1 1 _ l + (u/b){v + v- r ) + u 2 /b 2 

l-pL po (uv)l-pL po (u/v) ~ (1-p) 2 + (u/b)(v + v~ 1 )(l-p) + u 2 /b 2 

1 { V +(1- V(l -P))u 2 /b 2 

1-p (1 -p) 2 + (u/b)(v + v~ 1 )(l -p) +u 2 /b 2 ' 

This is nonnegative and completely monotone as a function of v + v _1 . 
Example 3.17. In the case of Example 3.12, the Levy measure of p, c is 

oo oo 
n— m— 1 

where a m = ^(q m — (—r/p)™). This Levy measure is not absolutely continuous. Thus, 
p c is never GGC for any parameters u, v, w and c. 



4. Exponential integrals for independent increment 
processes 

We say that a process X = {X t = X^=i t > 0} is a compound sum process, if the {X{\ 
arc i.i.d. random variables, {M t ,t> 0} is a renewal process and they arc independent. 
When {M t } is a Poisson process, X is nothing but a compound Poisson process and is 
a Levy process. Unless {M t } is a Poisson process, X is no Levy process. In this section, 
we consider the case when either £ or rj is a compound sum process and the other is an 
arbitrary Levy process. Although (£, rj) is not a Levy process, the exponential integral 
(1.1) can be defined and its distribution can be infinitely divisible and/or GGC in many 
cases as we will show in the following. 

Case 1: The process £ is a compound sum process 

First, we give a condition for the convergence of the exponential integral (1.1) when (£, rf) 
is not a Levy process. 
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Proposition 4.1. Suppose that (£t,?7i)t>o is a stochastic process where £ and rj are 
independent, r\ is a Levy process and £ t = * s a compound sum process with 

i.i.d. jump heights Xi,i = 1,2, ... , and i.i.d. waiting times Wi- Then (1.1) converges in 
probability to a finite random variable if and only if 

&->oo a.s. and [ ( -ijgg JPflq |edg)<oo (4.1) 
for Az{x)=J {0 x) P(X l >u)du. 

Proof. As argued in the proof of Proposition 3.2, we can rewrite the exponential integral 
as perpetuity 

[ e-C«-d^=f;(n^]S i , 
■ / (°.«0 £o\£=i / 

where A4 = c~ Xi and = r\ w . . By Theorem 2.1 of [6] the above converges a.s. to a finite 
random variable if and only if Y[7=i A% ~> a.s. and 

./(l.oo) V^( lo g?)/ 

for = Jl Q ^ P(— logAi > u)du. Using the given expressions for A\ and B\ in our 
setting, we observe that this is equivalent to (4.1). It remains to show that a.s. convergence 
of the perpetuity implies convergence in probability of (1.1). Therefore, remark that 

/ e~ is - dr) a = / c~ Cs - drj s + c~ Ctj »* (r? t - rj T ), 

J(0,t] J(P,T Mt ] 

where the first term converges to a finite random variable while the second converges in 
probability to since sup te[TM(iTVt+i) \r) t ~ Vr Mt I = sup te[0im) |?? t |. □ 

Now we can extend Proposition 3.2 in this new setting as follows. 

Proposition 4.2. Suppose that the processes £ and rj are independent and that = 
*l Xj is a compound sum process with i.i.d. jump heights Xi,i = 1,2,..., and i.i.d. 
waiting times W%,i = 1,2,..., such that (4-1) is fulfilled. Suppose that C(e~ Xl ) G H(M. + ) 
and £(r/ T ) G T(R + ) for r being a random variable with the same distribution as W\ and 
independent of rj. Then 

c( f e- £ '-dJer(t + ). 
V j (0,00) / 

Furthermore, if C(c~ Xl ) G H(WL), C(r) T ) G T(R) and C(ij T ) is symmetric, then C(V) G 
T(R). 
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In the following, we give some examples fulfilling the assumptions of Proposition 4.2. 

Example 4-3 (The case when r\ is nonrandom and C(W\) is GGC). For the 

case r)t = t, C{rj T ) belongs to T(R + ) if and only if C(t) does. Hence, for all waiting times 
which are GGCs and for a suitable jump heights of £, we have C(V) £ T(R+). 

Example 4-4 (The case when 77 is a stable subordinator and C(Wi) is GGC). 
Consider r\ to be a stable subordinator having Laplace transform L^(u) = cxp{— u"} 
with < a < 1. Then the Laplace transform of B := r\ T is given by Lb (it) = L T (ii Q ). This 
function is HCM if and only if r is GGC, since by Proposition 2.1, L r is HCM and hence 
also its composition with x a . Thus, whenever C(t) = C(Wi) is GGC, C(rj T ) is GGC, too, 
fulfilling the assumption of Proposition 4.2. 

Example 4-5 (The case when 77 is a standard Brownian motion and C(Wi) 
is GGC). Given that 77 is a standard Brownian motion, 771 has characteristic function 
E e i*vi =cxp(-z 2 /2), which yields Lb (it) = L T (it 2 /2). We can not see C(B) £ T(R+) 
from this, and in fact C(B) is in T(R) but not in T(R + ) (see Bondesson [4], page 117). 
Then using that 77 is symmetric, we can apply (4) in Proposition 2.1 and conclude that 
C{V) £ T(R) for suitable jump heights of \. 

Example 4-6 (The case when 77 is a Levy subordinator and C(W\) is a half 
normal distribution). The 1/2-stable subordinator 77 and the standard half normal 
random variable t have densities, respectively, given by, 

2a/7T V 7T 

These yield the density function of i] T as 

fn-r( x ) = J fvv ( x )fr (y) dy = ^ n 1 + x - 

Interestingly, this is an F distribution (see Sato [15], page 46) and since a random variable 
with an F distribution is constricted to be the quotient of two independent gamma 
random variables, we have that C(r) T ) £ T(R+). 

Case 2: The process t] t is a compound sum process 

Again we start with a condition for the convergence of (1.1). It can be shown similar to 
Proposition 4.1. 

Proposition 4.7. Suppose £ to be a Levy process and that rj t = X)i=i * s a compound 
sum process with i.i.d. jump heights Yi,i = 1,2,..., and i.i.d. waiting times Ui, i = 1, 2, . . . . 
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Then (1.1) converges a.s. to a finite random variable if and only if 



6 -> oo 



a.s. and 




( 



A/(log<7) 



log 9 



) 



P(|Yi| edg) <oo 



(4.2) 



for A n (x) = J {0 x) Pitu, > u) du. 



In the same manner as before, we can now extend Proposition 3.11 to the new setting 
and obtain the following result. 

Proposition 4.8. Let £ and r\ be independent and assume r\i = Xa=i ^ ^° be a compound 
renewal process with i.i.d. jump heights Yi,i = 1,2, ... , and i.i.d. waiting times Ui such 
that (4-2) holds. Suppose that C(Y\) G T(R + ) and C(e~^ T ) G H(M.+ ) for a random variable 
t having the same distribution as U\ and being independent o/£. Then 



Furthermore, if £(e"^) G H(M), £(Yi) G T(M) and £(Yi) is symmetric, then C(V) G 
T(R). 



The following is a very simple example fulfilling the assumptions in Proposition 4.8. 

Example 4.9 (The case when r\ is a random walk and £ is a standard Brownian 
motion with drift). Suppose £t — B t ■ + at is a standard Brownian motion with drift 
a > and U\ is degenerated at 1. Then £(e~^ T ) = C(e~ a e~ Bl ) is a scaled log-normal 
distribution and hence in H(M.+). So for all GGC jump heights C(Y\), the exponential 
integral is GGC. 
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